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ABSTRACT 

This paper, the 6'^ in the Semi-Empirical Analysis of Galaxies series, studies the evo- 
lution of 82302 star-forming (SF) galaxies from the Sloan Digital Sky Survey. Star 
formation histories (SFH) are derived from detailed spectral fits obtained with our 
publicly available spectral synthesis code starlight. Our main goals are to explore 
new ways to derive SFHs from the synthesis results and apply them to investigate how 
SFHs vary as a function of nebular metallicity {Zneb)- A number of refinements over 
our previous work are introduced, including (1) an improved selection criterion; (2) a 
careful examination of systematic residuals around H/3; (3) self-consistent determina- 
tions of nebular extinctions and metallicities; (4) tests with several Zneb estimators; 
(5) a study of the effects of the reddening law adopted and of the relation between 
nebular and stellar extinctions and the interstellar component of the Nal D doublet. 

Our main achievements may be summarized as follows: (1) A conventional corre- 
lation analysis is performed to study how global properties relate to Zneb, leading to 
the confirmation of previously known relations, such as those between Zneb and galaxy 
luminosity, mass, dust content, mean stellar metallicity and mean stellar age. (2) A 
simple formalism which compresses the results of the synthesis while at the same time 
yielding time dependent star formation rates (SFR) and mass assembly histories is 
presented. (3) A comparison of the current SFR derived from the population synthe- 
sis with that obtained from Ha shows that these independent estimators agree very 
well, with a scatter of a factor of two. An important corollary of this finding is that we 
now have a way to estimate SFR in galaxies hosting AGN, where the Ha method can- 
not be applied. (4) Fully time dependent SFHs were derived for all galaxies, and then 
averaged over six Zneb bins spanning the entire SF wing in the [OIIIJ/H/3 x [Nil] /Ha 
diagram. (5) We find that SFHs vary systematically along the SF sequence. Though 
all star-forming galaxies formed the bulk of their stellar mass over 1 Gyr ago, low 
Zneb systems evolve at a slower pace and are currently forming stars at a much higher 
relative rate. Galaxies at the tip of the SF wing have current specific SFRs about 
2 orders of magnitude larger than the metal rich galaxies at its bottom. (6) At any 
given time, the distribution of specific SFRs for galaxies within a Z„eb-bin is broad 
and approximately log-normal. (7) The whole study was repeated grouping galaxies 
within bins of stellar mass and surface mass density, both of which are more fun- 
damental drivers of SFH. Given the existence of strong Zneb~M*-^* relations, the 
overall picture described above remains valid. Thus, low (low E^,) systems are the 
ones which evolves slower, with current specific SFRs much larger than more massive 
(dense) galaxies. (8) This overall pattern of SFHs as a function of Zneb, or is 
robust against changes in selection criteria, choice of evolutionary synthesis models 
for the spectral fits, and differential extinction effects. 



Key Vifords: galaxies: evolution - galaxies: statistics - galaxies: stellar content. 
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1 INTRODUCTION 



The Sloan Digital Sky Survey (SPSS. IVork et alJ [ioooh . 
with its homogeneous spectroscopic and photometric data 
on hundreds of thousands of galaxies has revolutionized our 
perception of the world of galaxies in the local Universe. The 
enormous amount of objects allowed one to reveal trends 
that had not been suspected before. For example, while is 
was k nown since the work of [Baldwin. Phillips fc TerlevichI 
that objects ionized by massive stars and active 
galactic nuclei (AGN) live in different zones of emission line 
ratios diagrams, the fact that emis sion line galaxies are dis - 
tributed in two well defined wings l|Kauffmann et al.ll20033 ) 
in the famous [OIII]5007/H/3 vs [NII]6583/Hq; diagnostic di- 
agram (heareafter, the BPT diagram) came as a surprise. 

The left wing of the BPT diagram can be understood as 
a sequence in metallicity of normal star-forming (SF) galax- 
ies. The present-day nebular metallicity (Zneb) of a galaxy 
is intimately connected with its past star formation history 
(SFH). The main goal of this paper is to explore this link. 
The focus of many of the pioneering studies of SFH of SF 
galaxies was instead the variation of SFH with the Hub- 
ble type. ISearle. Sargent fc Bagnuolcl (119731 ) . for instance, 
assumed a simple model for the SFH and calculated UBV 
colors for simulated galaxies. Comparing simulated and ob- 
served colors, they concluded that morphological type alone 
does not explain the differences in SFH, proposing that the 
galaxy mass should also be used as a tracer of star forma- 
tion. 

[Gallagher HI. Hunter fc Tutukovl l| 19841 ) introduced a 
way to study the star formation rates (SFR) in three dif- 
ferent epochs of a galaxy's history. In order to achieve such 
time resolution, manifold indices were used: HI observations, 
dynam ical mass e s, _B- band and Ha luminosities, and UBV 
colors. ISaiidagel ll 19861) applied so me of the techniques pre- 
sented by Gallagher et al.l |l983) to investigate differences 
in SFH along the the Hubble sequence. I n the same vein, 
iKennicutt Jr.. Tamblvn fc CongdonI (|l994l ) derived the SFR 
for SF objects from the Ha luminosity and UBV colors, and 
found that the SFH differences for galaxies of the same Hub- 
ble type have a much stronger relation with their disk than 
with their bulge. iGavazzi et al. 112002,) measured the present 
and past SFRs of late-type galaxies in nearby clusters from 
Ha imaging and near infrared observations, and also derived 
the global gas content from HI and CO observations. Most 
of these studies had to rely on many different indices and 
observations in order to measure an instantaneous SFR, or 
at most a 2-3 age resolution star formation history for SF 
galaxies. 

Bical i 19881 ) introduced a method to reconstruct SFHs 
in greater detail by mixing the properties of a base of star 
clusters of various ages (t*) and metallicities (Z-,,). In its 
original implementation, this method used a set of 5-8 ab- 
sorption line equivalent widths as observables, a grid of clus- 
ters arranged in 35 combinations of and Z*, and a sim- 
ple parameter space exploration technique limited to paths 
through the plane constra ined by chemical evolution 

arguments (see also ISch midt ct al.l [l99ll : iBica et aL 1 11994 
ICid Fernandes et al.ll200ll ). Its application to nuclear spec- 
tra of nearby galaxies of different types revealed systematic 
variations of the SFH along the Hubble sequence. 

For over a decade the most attractive feature of Bica's 



method was its use of observed cluster properties, empir- 
ically bypassing the limitations of evolutionary synthesis 
models, which until recently predicted the evolution of stel- 
lar systems at spectral resolutions much lower than the data. 
This is no longer a problem. Medium and high spectral res- 
olution stellar libraries, as well as updates in evolutionary 
tracks have been incorporated into evolutionary synthesis 
models in the past few years. The current status of these 
models and their ingredients is amply discussed in the pro- 
ceedi ngs of the lAU Symposium 241 l|Vazdekis fc Peletieij 
[2003). 

These advances spurred the development of SFH re- 
covery methods which combine the non-parametric mixture 
approach of empirical population synthesis with the ambi- 
tious goal of fitting galaxy spectra on a pixel-by-pixel basis 
using a base constructed with these new generation of evolu- 
tionary synthesis models. Methods based on spectral indices 
have also benefitted from these new mod els, and produced 
an im press ive collection of resu l ts (e.g., Kauffmann et all 
I2003allbll3 : ICallazzi etall l2005l : iBrinchmann et all \2Qm . 



However, current implementations of these methods do not 
reconstruct detailed SFHs, although they do constrain it, 
providing estimates of properties such as mass-to-light ra- 
tios, mean stellar age, fraction of mass formed in recent 
bursts, and ratio of present to past SFR. 

The first SFHs derived from full spectral fits of SPSS 
galax ie s wer e carried out with the MQPEP (iPanter et al.l 
|2003| . l2007l : iMathis. Chariot fc BrinchmannI |2006| ) and 
STARLIGHT codes. MOPEP results have been recently re- 
viewed by Panter (2007), so we just give a summary of the 
results achieved with starlight. 

STARLIGHT itself was the main topic of the first paper 
in our Semi- Empirical Analysis of Galax ies series (SEAGal). 
In SEAgal I (ICid Fernandes et aLll2005l ) we have thoroughly 
evaluated the method by means of simulations, astrophys- 
ical consistency tests and comparisons with t he results ob- 
tain ed by independent groups. In SEAGal II l|Mateus et al.l 
I2OOQ) we have revisited the bimodality of the galaxy popu- 
latio n in terms of spectral synthesis products. In SEAGal III 
( Stasiriska et al. 2006), we combined the emission lines dug 
out and measured from the residual spectrum obtained af- 
ter subtraction of the synthetic spectrum with photoioniza- 
tion models to refine the criteria to distinguis h between nor- 
mal SF galaxies and AGN hosts. SEAGal IV l|Mateus et all 
12007. ) deals with environment effects, studied in terms of the 
relations between mean age, current SFR, density, luminos- 
ity and mass. 

Only in SEAGal V (|Cid Fernandes et al.l |2007| ) we 
turned our attention to the detailed time dependent infor- 
mation provide d by the synthesis. We have used the entire 
Pata Release 5 l|Adelman-McCarthv et al]|2007l ) to extract 
the population of SF galaxies, and study their chemical en- 
richment and mass assembly histories. It was shown that 
there is a continuity in the evolution properties of galax- 
ies according to their present properties: Massive galaxies 
formed most of their stars very early and quickly reached 
the high stellar metallicities they have today, whereas low 
mass (metal poor) galaxies evolve slower. These findings 
are in agreement with recent studies of the mass assem- 
bly of larges samples of g alaxies through the f ossil record of 
their stellar populations l|Heavens et al.l [20041 ). and of stud- 
ies of the distribution in small samples of galaxies (e.g.. 
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ISkillman et alj |2003| and references therein) , but the gen- 
erahty of the resuh apphed to the entire population of SF 
galaxies was shown for the first time. 

In the present paper, we aim at a more complete view 
of the properties of SF galaxies, and their variations along 
the SF sequence in the BPT diagram, improving and ex- 
panding upon the results only briefly sketched in SEAGal 
V. In particular, we discuss in depth time averaged values 
of quantities such as the SFR and the SFR per unit mass, 
as well as their explicit time dependence. 

The paper is organized as follows. Section 2 describes 
our parent sample and explain our criteria to define normal 
SF galaxies. This section also explains how we deal with 
extinction and how we estimate the nebular metallicity. In 
Section 3, we discuss the global properties of galaxies along 
the SF sequence in the BPT diagram. In Section 4, we then 
proceed to explain our formalism to uncover the explicit time 
dependence of such quantities as the star formation rate. In 
Section 5, we show that the current star formation rate as 
estimated by the most commonly indicator - the Ha lumi- 
nosity - compares with that obtained from our stellar pop- 
ulation synthesis analysis. In Section 6, we analyse the SFH 
along the SF sequence, binning galaxies in terms of their 
present-day nebular metallicity. We show that, despite the 
important scatter at any there is a clear tendency for 

the SFH as a function of Z^eb, in that in the most metal-rich 
galaxies most of the stellar mass assembly occurred very fast 
and early on, while metal poor systems are currently form- 
ing stars at much higher relative rates. We also compute 
mean SFHs binning the galaxies with respect to the stellar 
mass and the surface mass density, which are expected to 
better express the causes of the evolution of galaxies. Sec- 
tion [7] discusses possible selection effects and other caveats. 
Finally, Section |8] summarizes our main results. 



2 DATA 

The data analysed in this work was extracted from the SDSS 
Data Release 5 (DR5: lAdelman-McCarthv et al.|[2007l ). This 
release contains data for 582471 objects spectroscopically 
classified as galaxies, from which we have found ~ 1.6 per 
cent of duplicates, that is, objects with multiple spectro- 
scopic information in the parent galaxy catalog. 

From the remaining 573141 objects we have selected 
our parent sample adopting the following selection crite- 
ria: 14.5 < rur < 17.77 and z > 0.002. The magnitude 
range comes from the definition of the Main Galeixy Sam- 
ple, whereas the lower redshift limit is used to avoid inclu- 
sion of intragalactic sources. The resulting sample contains 
476931 galaxies, which corresponds to about 82 per cent of 
all galaxies with spectroscopic data gathered by SDSS and 
publicly available in the DR5. These limits imply a reduction 
by ~ 17% in the sample studied in SEAGal V. 



2.1 STARLIGHT fits 

After correcting for Galactic extinction (with the maps of 
ISchlegel. Finkbeiner fc David Il998l and the reddening law 
of ICardelh. Clavton fc MathisI 11989 . using Rv = 3.1), the 
spectra were shifted to the rest-frame, resampled to AA = 1 
A between 3400 and 8900 A, and processed through the 



STARLIGHT spectral synthesis code described in SEAGal I 
and II. 

STARLIGHT decomposes an observed spectrum in terms 
of a sum of simple stellar populations (SSPs), each of which 
contributes a fraction Xj to the flux at a chosen normal- 
ization wavelength (Ao = 4020 A). As in SEAGal II, we 
use a base of A^<- = 150 SSPs extracted from the mod- 
els of Bruzual fc CharlotI l|2003l . BC03), computed for a 
IChabrier (2003. ) initi al mass function (IMF), "Padova 1994" 
evolutionarv tracks ( ^ongi et al.lll993l: iBressan et al ] |l993l : 



iFagotto et alJll994a|[bl:lGirardi et al.lll996l '). and STELIB li- 
brarv ijLe Borgne et all l2003l l . The base components com- 
prise 25 ages between = 1 Myr and 18 Gyr, and 6 metal- 
licities, from Z*j = 0.005 to 2.5 solar. Bad pixels, emission 
lines and the NaD doublet are masked and left out of the 
fits. The emission line masks were constructed in a galaxy- 
by-galaxy basis, f ollowi ng the methodology outlined in SEA- 
Gal II and lAsaril l)2006l 'l. starlight outputs several physical 
properties, such as the present-day stellar mass, stellar ex- 
tinction, mean stellar ages, mean metallicities as well as full 
time dependent star formation and chemical evolution histo- 
ries, which will be used in our analysis. Section [4. 1 [ describes 
aspects of the code relevant to this work. 

Fitting half a million spectra represented a massive 
computational effort, carried out in a network of over 100 
computers spread over 3 continents and controlled by a 
specially designed PHP code. This huge d'atabase of spec- 
tral fits and related products, as well as starlight itself, 
are publicly available in a Vir tual Observatory environm ent 
at www.starlight.ufsc.br (see ICid Fernandes et al] |2007| . in 
prep.). 

Examples of the spectral fits obtained for 5 star-forming 
galaxies are shown in Fig [T] We have ordered the galaxies 
according to their nebular metallicity {Zneb), as defined in 
Section [2]4l to illustrate how spectral characteristics change 
along the Z„eb sequence. Metal-poor galaxies (top) show 
blue spectra and strong emission lines in comparison to the 
redder spectra and weaker emission lines of galaxies with a 
metal-rich ISM (bottom). The middle panels in Fig [T] show 
the fractional contribution to the total fiux at Ao — 4020 
A of simple stellar population (SSP) of age t*. These pan- 
els show that young stellar populations make a dominant 
contribution in galaxies with low Zneb, whereas at higher 
nebular metallicities a richer blend of stellar ages is present. 



2.2 Emission lines 

2.2.1 General procedure 

Emission lines were measured fitting gaussians to the resid- 
ual spectra obtained after subtraction of the stellar light 
using an updated version of the line-fitting code described 
in SEAGal III. The main transitions used in this study are 
up, [OIII]A5007, Ha and [NII]A6584. In the next section 
these lines are used to define the sub-sample of star-forming 
galaxies which will be studied in this paper. 



2.2.2 The special case of Hf3 

We find that a zero level residual continuum is adequate to 
fit the emission lines, except for H/3. Inspection of the spec- 
tral fits shows that the synthetic spectrum is often overesti- 
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Figure 1. Five examples of the spectral fits. Left panels show the observed (black) and fitted (red) spectra, both normalized at Aq = 4020 
A. Magenta lines mark regions not used in the fits either because they contain emission lines or because of artifacts in the data. Middle 
panels illustrate the fraction of light at Aq associated to each of the 25 SSP ages used in the fits. Curves represent a 0.5 dex smoothed 
version of the population vector. Right panels show SDSS 25.6" X 25.6" images (~ 12 X 12-34 X 34 kpc^). Galaxies in this plot are 
ordered according to their nebular metallicity {Zj^^i,; see Section l2.4l l. From top to bottom, Z„^i, = 0.29, 0.43, 0.61, 0.84 and 0.97 Zq. 



mated in the continuum around H/3, creating a broad, ~ 200 
A wide absorption trough in the residual spectrum. This 
problem, which can hardly be noticed in Fig[l] becomes ev- 
iden t when averaging m any residual spectra (see SEAGal V 
and IPanter et al-lTioOTl ) , and tends to be more pronounced 
for older objects. The comparison b etween STELIB stars 
and theoretical models presented by iMartins et al.l (|2005l ) 
gives a clue to the origin of this problem. A close inspection 
of Fig 21 in their paper shows that the STELIB spectrum 
has an excess of flux on both sides of H/3 when compared to 
the model spectrum. This suggests that the "H/3 trough" is 
related to calibration issues in the STELIB library in this 
spectral range. This was confirmed by starlight experi- 
ments which s howed that the problem disapp ears using the 
SSP spectra of iGonzalez-Delgado et ahl (|2005l . based on the 



Martins et al. 2 005 library) or those constru cted with the 
MILES hbrary l|Sanchez-Blazquez et aDl2006l )FI 

Though this is a low amplitude mismatch (equivalent 
width ~ 3 A spread over ~ 200 A), it makes H/3 sit in 
a region of negative residual flux, so assuming a zero level 
continuum when fltting a gaussian may chop the base of the 
emission line, leading to an underestimation of its flux. To 
evaluate the magnitude of this effect we have repeated the 
H/3 fits, this time adjusting the continuum level from two 
side bands (4770-4830 and 4890-4910 A). On average, the 
new flux measurement are 2% larger than the ones with the 
continuum flxed at zero. The difference increases to typically 
4% for objects with Whi3 < 5 A, but S/N > 3 in the line, 

^ We thank Drs. Enrique Perez, Miguel Cervino and Rosa 
Gonzalez-Delgado for valuable help on this issue. 
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Figure 2. The SF sample in the BPT plane, chopped into bins of 



nebular abundance Z^eb ■ 



(O/H) 



All lines have been corrected 



(0/H)q ■ 

by reddening (see Section [2.5l l. These same bins are used through- 
out this paper. The number of galaxies in each bin is given on the 
left. On the right, we show the corresponding mean Z^^i, and log 
mean values (in solar units). Galaxies close to bin borders are 
not plotted for clarity. 



and 7% for Whi3 < 2 A. Noise in the side bands introduces 
uncertainties in the measurement of the fiux, but at least it 
removes the systematic effect described above, so the new 
measurements should be considered as more accurate on av- 
erage. We adopt these new H/3 measurements throughout 
this work. Using the zero continuum measurements changes 
some of the quantitative results reported in this paper min- 
imally, with no impact on our general conclusions. 



2.3 Definition of the Star Forming Sample 

Since the pioneering work of iBaldwin. PhiUips fc TerlevichI 
l|l98ll ). emission line objects are classified in terms of 
their location in diagrams involving pairs of line ra- 
tios. As explained in SEAGal III, the [NII]A6584/Ha vs. 
[OIII]A5007/H/3 diagram (the BPT diagram) is the most 
useful for this purpose, mainly due to the partially secondary 
nature of N (i.e., t he increase of N/Q as O/H increases, e.g., 
iLiang et al.ll2006l : iMolla et al.ll2006l ). 

Our sample of star-forming galaxies is composed by 
objects which are below the line separating normal star- 
forming galaxies and AGN hosts proposed in SEAGal III. 
We have imposed a lower limit of 3 in S/N on the 4 lines 
in the BPT diagram, and S/N > 10 in the 4730-4780 A 
continuum to constitute our main sample (the SF sample). 
The 82302 galaxies composing the SF sample are shown in 
Fig. [Hon the BPT plane. 

Both the STARLIGHT fits (and thus aU SFH-related pa- 
rameters) and the emission line data are afTected by the 
quality of the spectra. To monitor this effect, we have de- 
fined a "high-quality" sub-set (the SF*^' sample) by doubling 
the S/N requirements for the SF sample, i.e., S/N > 6 in 
all 4 lines in the BPT diagram and a continuum S/N of 20 
or better. A total of 17142 sources satisfy these criteria. 

Fig.[3]shows the distributions of observational and phys- 



ical properties for the samples. Naturally, the SF'"' sample 
is skewed towards closer and brighter galaxies with respect 
to the SF sample, but in terms of physical properties such 
as stellar mass, mean age and nebular metallicity the two 
samples are similar. 



2.4 Nebular Metallicity Estimate 

It is well known that the SF-wing in the BPT diagram is a 
sequence in nebular metallicity (SEAGal III and references 
therein), which we quantify by the oxygen abundance ob- 
tained throu gh the O3N2 = [OIII]5007/[NII]6583 index as 
calibrated by IS tasihska l (j200d) : 



log Z„eb = log 



(O/H) 
(O/H), 



-0.14 -0.25 log O3N2 



(1) 



4.9 X 10" 



where we have adopted (0/H)„ = 
dAllende Prieto. Lambert fc Asplundl^2001^ . 

We have chosen to use the O3N2 indicator to estimate 
the average oxygen abundance in the ISM of SF galaxies 
mainly because it is a single-valued indicator and it can 
be easily related to the position of galaxies in the classi- 
cal BPT diagram. However, this indicator is affected by the 
the presence of diffuse ionized gas in galaxies and by the 
fact that the N/ O ratio depends on the way chemical ev olu- 
tion proceeded (|Chiappini. Romano fc Matteuccil |2003| ) . In 
addition, for the lowest metallicity galaxies, O3N2 is not 
sensitive to O/H anymore, as a wide range of metallicities 
corr espond to t he sam e value of O3N2, as can be seen in Fig. 
3 of IStasihskal (|2006l '). From that figure, the O/H given by 
equation UJ lies towards the upper end of the possible val- 
ues of O/H. We considered using the [Arlll] 7135 / [OII I] 5007 
(Ar303) index which, as argued by IStasiiiskal (|2006h . does 
not suffer from the problems mentioned for the O3N2 index. 
It turns out that, in objects where the [Arlll] line could be 
measured, ArsOs and O3N2 are extremely well correlated 
(with a Spearman correlation coefficient of Rs = 0.58). Un- 
fortunately, the quality of the SDSS spectra did not allow 
us to measure the [Arlll] line intensity with sufficient ac- 
curacy in a large number of objects, and principally in the 
zone where it would have b een helpful to break the O 3N2 vs 
O/H degeneracy. Using the lPilvugin fc ThuanI (|2005h metal- 
licity calibration based on [OIII]5007/H/3 and [OII]3727/H/? 
adds only a tiny fraction of galaxie s. The same applies when 
using the O/H values obtained by llzotov et al.l (I2OO6I ) from 
direct methods using the electron temperature derived from 
[OIII]4363/[OIII]5007. When comparing our measures with 
theirs, there is a systematic offset of 0.2 dex and a rms of 
0.13 dex in Z„eb for the 177 objects we have in common. 
This effect is greater the lower Znet is. We thus decided to 
use O3N2 as a nebular metallicity indicator all along the 
SF galaxy sequence, keeping in mind that equation ((T| will 
tend to attribute a metallicity Zneb of about 0.2 Zq for the 
galaxies with the lowest observed O3N2 in our sample. 

Other metal licity estimates have b een used for galax- 
ies. For example, iTremonti et al.l (|2004l ) obtained the nebu- 
lar metallicities by comparing the observed line ratios with 
a large data base of photoionization models. While a priori 
ap pealing, th i s meth od is not devoid of problems, as shown 
by lYin et aU (|2007l '). There is a systematic offset of -0.28 
dex and a rms of 0.09 between our nebular metallicities and 
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Figure 3. Normalized histograms of observed and physical properties for the SF (solid), SF'"' (dotted) samples. 



theirs. Their method also yields a larger range of values 
for Zneb- For the SF sample, their calibration covers from 
Zneb ~ 0.78 to 2.70 Zq for the 5 to 95 percentile ranges, 
whereas our calibration covers from 0.47 to 1.13 Zq for the 

same percentile ranges. 

The calibration bv lPettini fc Pagell l|2004l ) IS more sim- 
ilar to our own. There is a slight offset of -0.04 dex with re- 
spect to our calibration and the dispersion for the SF sample 
is 0.03 dex. Their calibration also stretches a little the Zn^t 
ranges: Z„^t = 0.45 to 1.36 Zq for the 5 to 95 percentile 
ranges. 

Although we believe that our calibration is likely more 
reliable, we have performed all the computations in this pa- 
per al so with the Trcmoriti ct al. (2004) and Pcttini & Pagcl 
(|2004l ) calibrations. While the results differ in absolute 
scales, the qualitative conclusions remain identical. 



2.4-1 Zneb bins 

As seen above, both physical and mathematical motivations 
make O3N2 a convenient index to map galaxy positions 
along the SF wing in the BPT diagram. From equation ^ 
one sees that a given value of Zmb using this index defines 
a straight line of unit slope in the BPT diagram. 

Our SF sample spans the Zn^b = 0.2-1.6 Zq range from 
the tip of the SF-wing to its bottom. In Fig. [2] this interval 
is chopped into 6 bins of width Alog Z„eb ~ 0.13 dex, ex- 
cept for the one of lowest metallicity which is twice as wide 
to include more sources. Table [T] lists some properties of 
galaxies in each of these bins, which are hereafter labeled 
A-F. Galaxies inside these bins will be grouped together in 
the analysis of star-formation presented in Section |6] Note 
that the bias in the determination of Z^eb from O3N2 at low 
metallicities has no consequence for our study, since almost 
all the objects from bin A remain in this bin. 



2.5 Extinctions 

2.5.1 Stellar extinction 

As explained in SEAGal I, starlight also returns an esti- 
mate of the stellar visual extinction. Ay, modeled as due to 
a foreground dust scree n. This is obviously a simplification 
of a complex problem (|Witt. Thronson Jr. fc Capuano Jrl 
1 19921 ). so that Av should be called a dust attenua- 
tion parameter instead of extinction, although we do not 
make this distinction. Previous papers in this series have 



used the ICardelli. Clayton fc Mathi^ (|l989l , CCM) redden- 
ing law, with Rv = 3.1. In order to probe different 
recipes for dust attenuation, we have selected 1000 galax- 
ies at random from the SF'"' sample and fitted them 
with four other functions: the starbur s t atte nuation law of 
ICalzetti. Kinney fc Storchi-BergmannI (|l993 ). the SMC and 
LMC c urves from ICordon et al.| (|2003l ) and the A"'' '^ law 
used by lKauffmann et al.l ( 200331 ). 

We find that the quality of the spectral fits remains 
practically unchanged with any of these 5 laws. Averaging 
over all galaxies the SMC law yields slightly better x^'s, 
followed closely by the Calzetti, A"°-^, LMC, and CCM, in 
this order. As expected, these differences increase with the 
amount of dust, as measured by the derived Ay values or by 
Ha/H/3. Yet, KS-tests showed that in no case the distribu- 
tions of x^'s differ significantly. This implies that the choice 
of reddening law cannot be made on the basis of fit quality. 
A wider spectral coverage would be needed for a definitive 
empirical test. 

When using different recipes for the dust attenuation, 
the synthesis algorithm has to make up for the small vari- 
ations from one curve to another by changing the popula- 
tion vector and the value of Ay. To quantify these changes 
we compare results obtained with the Calzetti and CCM 
curves, and consider only the most extincted objects. Com- 
pared to the results for a CCM law, with the Calzetti law 
the mean stellar age decreases by 0.09 dex on the median 
(q ualitatively in agre ement with the results reported in Fig 6 
of lPanter et al.|[2007l ). the mean stellar metallicity increases 
by 0.05 dex. Ay increases by 0.07 mag and stellar masses 
decrease by 0.02 dex. These differences, which are already 
small, should be considered upper limits, since they are de- 
rived from the most extincted objects. Somewhat larger dif- 
ferences are found when using the SMC and A"" '^ laws. For 
instance, compared to the Calzetti law, the SMC law pro- 
duces mean stellar ages 0.15 dex younger and masses 0.07 
dex smaller, again for the most extincted objects. 

We have opted to use the Calzetti law in our starlight 
fits and emission line analysis. The reasons for this choice 
are twofold: (1) The Calzetti law yields physical properties 
intermediate between the SMC and CCM laws, and (2) this 
law was built up on the basis of integrated observations of 
SF galaxies, similar to the ones studied in this paper. In any 
case, the experiments reported above show that this choice 
has little impact upon the results. 

STARLIGHT also allows for population dependent extinc- 
tions. Although tests with these same 1000 galaxies sample 
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Bin A 


Bin B 


Bin C 


Bin D 


Bin E 


Bin F 


logZ^eb min [Z©] 


-0.710 


-0.450 


-0.320 


-0.190 


-0.060 


0.070 


log Zr,eb P05 


-0.608 


-0.435 


-0.311 


-0.180 


-0.053 


0.071 


log Zneh P50 


-0.494 


-0.364 


-0.242 


-0.111 


-0.001 


0.081 




-0.454 


-0.324 


-0.194 


-0.064 


0.054 


0.112 




-0.450 


-0.320 


-0.190 


-0.060 


0.070 


0.200 


logM^ p05 [Mq] 


7.146 


7.938 


8.693 


9.277 


9.800 


10.089 




8.319 


8.922 


9.452 


9.958 


10.460 


10.678 
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(logi*)z, p05 [yr] 


6.755 


7.428 


7.695 


7.915 


8.127 


8.167 


(logt*>i p50 


7.762 


8.155 


8.385 


8.567 


8.726 


8.724 


{logt*)i p95 


8.556 


8.900 


9.139 


9.255 


9.276 


9.255 


logl-Ha p05 [Lq] 


5.167 


5.139 


5.492 


5.988 


6.551 


6.956 


logLifa p50 


6.475 


6.475 


6.750 


7.149 


7.544 


7.816 


logLi/c« p95 


7.834 


7.736 


7.961 


8.172 


8.426 


8.559 


log b p05 


-0.157 


-0.392 


-0.698 


-0.947 


-0.973 


-0.890 


log b p50 


0.769 


0.436 


0.211 


-0.006 


-0.213 


-0.199 


log b p95 


1.527 


1.192 


0.905 


0.644 


0.366 


0.272 



Table 1. Statistics of properties in bins A— F. 



show that in general one obtains larger Ay for young popula- 
tions, as expected, simulations show that, as also expected, 
this more realistic modeling of dust effects is plagued by 
degeneracies which render the results unreliable (see also 
IPanter et all l2007l : SEAGal V). We therefore stick to our 
simpler but more robust single Ay model. 



2.5.2 Nebular extinction 

The nebular V-band ex tinction was compu ted from the 
Ha/H/3 ratio assuming a lCalzetti et al.l (|l994l ) law: 



= 7.96 log (Ha/H/3U. 



(Ha/H/3)i„ 



(2) 



where (Ha/H/3)obs and (Ha/H/3)int are the observed and 
intrinsic ratio respectively. Instead of assuming a constant 
value, we take into account the metallicity dependence of 
(HQ/H/3)int, which varies between 2.80 and 2.99 for Z^eh in 
the 0.1 to 2.5 Zq range, as found from the photoionization 
models in SEAGal Ilfl 

We obtain the intrinsic ratio as follows. We start by 
assuming (Ha/H/3)int = 2.86, from which we derive a first 
guess for A^-^'. We then use the dereddened [OIII]A5007 and 
[NII]A6584 line fluxes to calculate Zneb (eq. [l]). From our 
sequence of photoionization models (SEAGal III) and Zneb, 
we derive a new estimate for (HQ:/H/3)int, and hence Ay''^ 
(eq. [21 . It takes a few iterations (typically 2-3) to converge. 

For 1.6% of the objects (1.1% for the SF'" sample), 
Ha/H/3 is smaller than the intrinsic value, which leads to 
A"^"*" < 0. In such cases, we assume A^^ = 0. We have 
corrected both the [0III]/H/3 and [Nil] /Ha line ratios for 
dust attenuation for the remainder of our analysis. 



^ Note that the models take into account collisional excitation of 
Balmer lines, so that at low metallicities the intrinsic Ha/H/3 is 
different from the pure case B recombination value. 
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Figure 4. (a) Relation between stellar and nebular extinctions 
for the SF sample. Lines indicate the 5, 50 and 95 percentiles, (b) 
Equivalent width of the ISM component of the NaD doublet (as 
measured from residual spectra) against the stellar extinction. 



We find that Ay'' and Ay axe strongly correlated, as 
shown in Fig 3^. A robust linear fit including all points 
yields A^'^'' = 0.34 + 2.28^^. The ionized gas thus suffers ~ 
twice as much extinction as t he stellar continuum, c orrob- 
orating the results reported in lStasihska et aP l)2004h (with 
a different methodology) and SEAGal I (obtained with a 
smaller sample, different version of starlight and a CCM 
extinction curve), and in agreement with d etailed studies 
of nearby SF-galaxies (ICalzetti et al.|[l994l ). We also find 
that the difference between nebular and stellar extinctions 
increases systematically as the mean age of the stellar pop- 
ulation increases. 

Given the spatial association of the line emitting gas 
and the ionizing populations, these results ultimately im- 
ply a breakdown of our simple single-A^ modelling. In 
fact, STARLIGHT experiments with population dependent Ay 
point in the same direction, i.e, the need to allow young pop- 
ulations to suffer more extinction than older ones. To eval- 
uate to which extent this simplification affects the results 
reported in this paper. Sec. 17.31 presents experiments where 



8 Asari etal. 



the extinction of < 10 Myr components is set according 
to the empirical relation Ay^''(Ay) found above. 

2.5.3 Interstellar absorption as traced by the NaD doublet 

The most conspicuous spectroscopic feature of the cold ISM 
in the optical range is the NaD doublet at AA5890,5896 A. 
For a constant gas to dust ratio, the strength of this feature, 
which measures the amount of cold gas in front of the stars, 
should correlat e with A^, as found for far-lR bright star- 
burst galaxies (|Heckman et al.l l200(j). To perform this test 
for our sample, we first measure the flux of the NaD doublet 
in the residual spectrum, integrating from 5883 to 5903 A. 
We thus remove the stellar component of this feature, which 
is also present in stellar atmospheres, pa r ticularly late typ e 
stars (jjacobv. Hunter fc ChristianI Il984l : iBica et alj Il99ll ). 
In principle, this is a more precise procedure than estimat- 
ing the ste llar NaD from its rela t ion to other stellar absorp - 
tion hues l|Heckman et al. I I2OOOI : ISchwartz fc Martini [2OM ) ■ 
but since the NaD window was masked in all fits (precisely 
because of its possible contamination by ISM absorption), 
the stellar NaD predicted by the fits rely entirely on other 
wavelengths, so in practice this is also an approximate cor- 
rection. The residual flux is then divided by the continuum 
in this range (deflned as the median synthetic flux in the 
5800-5880 plus 5906-5986 A windows), yielding the excess 
equivalent width AWnbD, which says how much stronger 
(more negative) the NaD feature is in the data with respect 
to the models. 

FigHb shows the relation between AWNaD and Ay. The 
plot shows that these two independently derived quantities 
correlate strongly. Intriguingly, AWNaD converges to 0.8 A 
in the median as Ay — > 0. We interpret this offset from 
AWNaD = as due to the fact that the stars in the STELIB 
Hbrary have a Galactic ISM component in their NaD lines. 
This propagates to our spectral models, leading to an over- 
prediction of the stellar NaD strength, and thus AWNaD > 
when the ISM absorption approaches zero. 

Regardless of such details, the discovery of this astro- 
physically expected correlation strengthens the confidence in 
our analysis. Furthermore, it opens the interesting prospect 
of measuring the gas-to-dust ratio and study its relation with 
all other galaxy properties at hand, from nebular metallic- 
ities to SFHs. This goes beyond the scope of the present 
paper, so we defer a detailed analysis to a future communi- 
cation. 



3 CORRELATIONS WITH NEBULAR 
METALLICITY 

Galaxy properties change substantially from the tip of the 
SF-wing, where small, metal-poor Hll-galaxy-like objects 
live, to its bottom, populated by massive, lumino us galaxies 
with large bulge-to-disk ratios and rich in metals (jKennicuttl 
1 19981 ). The simplest way to investigate these systematic 
trends is to correl ate various p ropertie s with the nebu- 
lar m etallicity (e.g:.. lTremonti et a l. 2004; Brinchmann et al] 
I2OO4I ). 

In this section we correlate Z^eb with both observed 
and physical properties extracted from our stellar popula- 
tion fits. This traditional analysis, based on current or time- 



averaged properties, helps the interpretation of the more de- 
tailed study of time-dependent SFHs presented in the next 
sections. In fact, this is the single purpose of this section. 
Since most of the results reported in this section are already 
know or indirectly deducible from previous work, we will 
just skim through these correlations. 

Fig [5^ shows Z„eb against absolute r-band magnitude. 
This is the luminosity-nebular metallicity relation, previ- 
ously studied by many authors, and interpreted in terms 
of a mass-metallicity relation. Fig [5]d shows our version of 
the Mi,-Zrieb relation. Because of the expected bias in our 
Zncb estimate at the lowest metallicities (see Section |2.4|) . 
we expect the real mass-metallicity relation to be flatter at 
low M-t than seen i n this plot. 

As shown by iKauffmann et al.1 (|2003bl ). stellar mass 
and stellar surface mass density (E*) are very strongly re- 
lated. It is thus no surprise to find that E* also corre- 
lates with ZnfLb, as show n in Fig[5h. Our definition of is 
the same as adopted bv IKauffmann et al.l (l2003bl ). namely 
E* = Mt/27rr|o,z) where rso.z is the half light Petrosian 
radius in the z-band. 

Fig[5ji shows how nebular extinction increases system- 
atically with Zneb. One factor which surely contributes to 
this relation is the rise in dust grain formation with increas- 
ing g as metallicity, but oth er factors may come into play as 
weU (IStasiiiska et al.ll2004l '). 

Fig [5^ shows how nebular and stellar metallicities cor- 
relate. This important relation, first presented in SEAGal I 
(for different sample, SSP base and Z^eb scale), shows that 
stellar and ISM chemical enrichment levels scale with each 
other, as one would expect on the basis of simple chemical 
evolution scenarios. The large scatter in Fig [5^ is mostly in- 
trinsic (as we verified comparing the relation obtained for 
data of different qualities), in qualitative agreement with 
the idea that stellar and nebular metallicities reflect differ- 
ent evolutionary phases and react differently to the several 
processes which regulate the chem ical evolution of galax ies. 
A similar relation was obtained bv lGallazzi et al.l l|2005h us- 
ing different methods to estimate both stellar and nebular 
abundances. Even though we express both quantities in solar 
units, these two metallicities are derived by such radically 
different means that, as discussed in SEAGal V, they should 
not be compared in quantitative terms. 

The relation between the mean stellar age {logt*)L and 
Zneb, shown in Fig[Sf, reflects the fact that young stars have 
a larger share of the light output at the tip of the SF wing 
than at its bottom, where old populations have a greater 
weight. Metal-rich SF galaxies thus have a more continuous 
star-forming history than metal-poor ones, which are often 
dominated ( in light, but not m ass) by the latest generation 
of stars (e.g:.. ICorbin et al]|2006l ). This is another way to look 
at metallicity-age trend, discussed previously in the analysis 
of Fig. [1] Ultimately, this relation represents a summary of 
chemical evolution, in the sense that more evolved systems 
have a more enriched ISM. In a related vein. Fig [5^ shows 
how the ratio of current to mean past SFR (defined in Sec- 
tion l5.3p varies along the metallicity sequence of SF galaxies. 
This indicates that the lower-metallicity galaxies are slower 
in forming stars. When one considers the mass-metallicity 
relation (FiglSJj), this is just another way of looking at the 
downsizing effect (Heavens et al1l2004l : iThomas et al. I I2OO5I : 
iMateus et al.|[200i ). Finally, Fig[5]i shows the relation be- 
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Figure 5. Correlations of and (a) tlie absolute r-band magnitude, (b) the stellar mass, (c) the surface mass density, (d) the nebular 
extinction, (e) the mean stellar metallicity, (f ) the mean stellar age, (g) the ratio of current to mean past SFR, and (h) the Ho luminosity. 
Numbers in each panel report the Spearman rank correlation coefficient, and the lines mark the 5, 50 and 95% percentiles of 25 bins, 
3292 points in each bin. The right hand scale in panel (h) is also logLj^Q, but in units of 2 X 10* Lq, such that it can also be read as an 
estimate of the current SFR in M0yr~^ (Section [Sjl . 



tween reddening corrected Ha luminosity and Zne.b- The y- 
axis is given in units such that the values correspond ~ 
to the current SFR in M0 yr~^ (see Section [S}. The cor- 
relation, although statistically unquestionable, has a large 
scatter. This implies that galaxies in the 6 .Zneb bins defined 
in Fig [5] have heavily overlapping current SFRs. Section [S] 
presents independent confirmation of this fact. 

As expected, all correlations discussed above are also 
present for the SF'"' sub-sample. For most they are in fact 
somewhat stronger, whereas for samples defined with less 
stringent criteria the correlation strengths weaken, indicat- 
ing that noise in the data is responsible for part of the scat- 
ter in these relations. Finally, as is widely known and can 
be deduced from Fig [S] itself , there are many inter-relations 
between galaxy properties. Our use of Z„e6 as the "inde- 
pendent" variable axis in Figgis not meant to indicate that 
Zneb is the underlying cause of the correlations; it simply re- 
flects our interest in mapping physical properties of galaxies 
along the SF wing of the seagull in the BPT diagram. 



4 METHODS TO INVESTIGATE STAR 
FORMATION HISTORIES 

The main goal of this paper is to study how the SFH 
varies among SF galaxies. Most other investigations in 
this same line used absorption, emission or continuum 
spectral indices such as the 4000 A break, the H5 ab- 
sorption, the K, G and Mg bands, or the Ha luminos- 
ity and equiva l ent w i dths to cha r acter iz e the SFH (e.g., 
Raimann et al. ' I 2OOOI: iKong et al.l |2003|; iKauffmann et al.l 
2003b; Cid Fcrna ndes et all l2o"o3riBrinchmann et al.ll2004 
Westera et al.ll20oi r Our approach, instead, is to infer the 



SFH from detailed pixel-by-pixel fits to the full observed 
spectrum, thus incorporating all available information. 

Whereas our previous work concentrated on the first 
moments of the age and metallicity distributions, here we 
present some basic formalism towards a robust description 
of SFHs as a function of time. From the point of view 
of methodology, these may be regarded as "second-order" 
products. Astrophysically, however, recovering the SFH of 
galaxies is of prime importance. SEAGal V presented our 
first results in this direction, including empirically derived 
time-dependent mean stellar metallicities. In this section we 
expand upon these results, exploring new ways to handle 
the output of the synthesis, focusing of the SFHs. 



4.1 Compression methods 



As reviewed in Section 12.11 starlight decomposes an ob- 
served spectrum in terms of a sum of SSPs, estimating the 
Xj (j = 1 ■ ■ ■ Ni,) fractional contribution of each popula- 
tion to the flux at Aq — 4020 A. For this work we used 
a base of A^^ = 150 SSPs from BC03, spanning 25 ages 
between t*j — 1 Myr and 18 Gyr, and 6 metallicities 
(0.005 < Zi, j < 2.5Zq). Example fits were shown in Fig 

m 

Not surprisingly, the 150 components of the population 
vector (x) are highly degenerate due to noise, and astro- 
physical plus mathematical degene r acies, as confirmed by 
simulations in I Cid Fernandes et all (|2004l ') and SEAGal I. 
These same simulations, however, proved that compressed 
versions of the population vector are well recovered by the 
method. 

Different compression appr oaches exist among spectral 
synth e sis codes. In MOPED llHeavens. Jimenez fc Lavahl 
I2OOOI : iReichardt erall I2001I : fpanter et al.l |2003| . ^OOtT 
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iHeavens et al ] |2004l ). for instance, compression is done a pri- 
ori, replacing the full spectrum by a set of numbers associ- 
ated to each of the A'^* + 1 parameters (the mass fractions 
and metallicities i n several tim e bins plus a dust parame- 
ter). STECMAP (|Ocvirk et al.l liooe') performs a compres- 
sion by requiring the resulting SFH to be relatively smooth. 
The preference for a smooth solution over a ragged one is 
effectively a prior, but the algorithm adjusts the degree of 
smoothing in a data driven fashion, so we may call it an 
"on the fly" c ompression method . The same can be said 
about VESPA (|Toieiro et al.ll200'<1 ). a new code which com- 
bines elements from these two approaches, starlight is less 
sophisticated in this respect. Its only built-in compression 
scheme is that the final stages of the fit (after the Markov 
Chains reach convergence) are performed with a reduced 
base comprising the subset of the original A*"* populations 
which account for > 99% of the light. For our parent sample 
of 573141 galaxies the average size of this subset is N^^ = 24 
populations, while for the 82302 galaxies in the SF sample 
jycff _ (This difference happens because the full sample 
has many old, passive systems, which require relatively few 
SSPs, whereas SF galaxies have more continuous SF regimes, 
thus requiring more SSPs to be adequately fit.) Compres- 
sion beyond this level must be carried out a posteriori by 
the user. As explained in the next section, in this study we 
in fact compress this information into only four age bins by 
smoothing the population vectors. 

Previous papers in this series have taken this a posteri- 
ori compression approach to its limit, condensing the whole 
age distribution to a single number, the mean stellar age; 



(lo! 



■,U)l ^^Xj\ogt^,,j 



(3) 



where the subscript L denotes a light- weighted average. 
Mass-weighted averages are readily obtained replacing x 
by the mass-fraction vector jj. Similarly, stellar metallici- 
ties were only studied in terms of their mass- weighted mean 
value: 



(4) 



Simulations show that both of these quantities have 
small uncertainties and essentially no bias. Regarding prac- 
tical applications, these first moments proved useful in the 
study of several astrophysical relations, some of which have 
just been presented in Section |3] (see Fig. [5]). Notwithstand- 
ing their simplicity, robustness and usefulness, these aver- 
ages throw away all time dependent information contained 
in the population vector, thus hindering more detailed stud- 
ies of galaxy evolution. In what follows we explore novel 
ways to deal with the population vector which circumvent 
this limitation. 



4.2 Star Formation Rate as a Function of Time 

One alternative to characterize higher moments of the SFH 



is to b in x onto a ge-groups, a strategy t hat goes back tolBica 
119881 : see also ISchmidt et al. 1 Il99ll : [cid Fernandes et al 



define bin-limits, and produces a discontinuous description 
of the SFH. 

A method which circumvents these disadvantages is to 
work with a smoothed version of the population vector. We 
do this by applying a gaussian filter in log t*, with a FWHM 
of 1 dex. Given that our base spans ~ 4 orders of magnitude 
in t*, this heavy smoothing is equivalent to a description in 
terms of ~ 4 age groups, but with the advantage that Xs can 
be sampled continuously in logt*. This approach is analo- 
gous to smoothing a noisy high-resolution spectrum to one 
of lower resolution, but whose large-scale features (colours, 
in this analogy) are more robust. From the results in SEA- 
Gal I, where it was shown that 3 age groups are reliably 
recovered, we expect this smoothing strategy to be a robust 
one. Furthermore, averaging over a large number of objects 
minimizes the effects of uncertainties in the smoothed SFH 
for individual galaxies. 

Technically, the issue of age resolution in population 
synthesis is a complex one. lOcvirk et aLl(|2006l ). for instance, 
find that bursts must be separated by about 0.8 dex in log 
to be well distinguished from one another with their SFH 
inversion method. Considering that the age range spanned 
by our base is 4.2 dex wide, one obtains 5 "independent" 
time bins. This number is similar to t hat (6 bins) used by 
iMathis. Chariot fc BrinchmannI l|2006l ) to describe the SFH 
of SDSS gala xies (cover ing a wider A-range than those sim- 
ulated by Ocvirk et al.l but at lower S /N) with a variant of 
the MOPED code. Up to 12 time bin s were used in other ap- 
phcations of MOPED. IPanter et al] (jiooj) argue that this 
may be a little too ambitious for individual galaxies, but 
uncertainties in this overparameterize d description average 
out in applications to large samples. iToieiro et al.l (|2007f ) 
have a useful discussion on the number of parameters that 
can be recovered with synthesis methods. By using their 
VESPA code and calculating the number of parameters on 
the fly for each individual object, they find that tipically 2- 
8 parameters can be robustely recovered for SDSS spectra. 
Hence, despite the complexity of the issue, there seems to be 
some general consensus that the age resolution which can be 
achieved in practice is somewhere between 0.5 and 1 dex, so 
our choice of smoothing length is clearly on the conservative 
side. 

A further advantage of this continuous description of 
the SFH is that it allows a straight-forward derivation of 
a star- formation rate (SFR). Recall that we describe a 
galaxy's evolution in terms of a succession of instantaneous 
bursts, so a SFR is not technically definable unless one as- 
sociates a duration to each burst. The SFR(fi) function is 
constructed by sampling the smoothed mass-fraction vector 
jig in a quasi-continuous grid from logt* = 5.6 to 10.5 in 
steps of Alogt* = 0.1 dex, and doing 



SFR(t,) = 



dM;{t^) AM^{t^) M^loge fi^U) 



dt^. 



At-, 



A log t-t. 



(5) 



I2OOII ). Though useful, this approach introduces the need to 



where is the total mass converted to stars over the galaxy 
history until ti, — 0, and /is(t*) is the fraction of this mass 
in the bin0 

^ The superscript c is introduced to distinguish from the 
mass still locked inside stars (M*), which must be corrected for 
the mass returned to the ISM by stellar evolution. This distinc- 
tion was not necessary in previous SEAGal papers, which dealt 
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SSFR(t 



We can also define the time dependent specific SFR: 
1 dM^iU) _ loge fi'=,{U) 



MS 



U A log U 



(6) 



which measures the pace at which star-formation proceeds 
with respect to the mass already converted into stars. This is 
a better quantity to use when averaging the SFH over many 
objects, since it removes the absolute mass scale dependence 
of equation ([5]). 

Three clarifying remarks are in order. (1) All equa- 
tions above are marginalized over Zi,, i.e., SFR(f*) = 
y^^y SFR(f*, Zi,) measures the rate at which gas is turned 
into stars of any metallicity. (2) The upper age limit of our 
base (18 Gyr) is inconsistent with our adopted cosmology, 
which implies an 13.5 Gyr Universe. Given the uncertain- 
ties in stellar evolution, cosmology, observations and in the 
fits themselves, this is a merely formal inconsistency, and, 
in any case, components older than 13.5 Gyr can always be 
rebinned to a cosmologically consistent time grid if needed. 
(3) Finally, since our main goal is to compare the intrin- 
sic evolution of galaxies in difi'erent parts of the SF wing in 
the BPT diagram, throughout this paper we will consider 
ages and lookback times in the context of stellar-evolution 
alone. In other words we will not translate t* to a cosmolog- 
ical lookback time frame, which would require adjusting the 
scale by adding the z-dependent lookback time of each 
galaxy. 

4.3 Mass Assembly Histories 

Another way to look at the population vector is to compute 
the total mass converted into stars as a function of time: 



(7) 



which is a cumulative function that grows from to 1, start- 
ing at the largest f*, tracking what fraction of A/^ was con- 
verted to stars up to a given lookback time. 

We sample rjl in the same logt* = 5.6-10.5 grid used to 
describe the evolution of the SFR, but here we operate on 
the original population vector, not the smoothed one. Since 
most of the mass assembly happens at large t*, computing ryj 
with the smoothed SFHs leads to too much loss of resolution. 
In essence, however, »75(t*) and SSFR*(tt) convey the same 
physical information in different forms. 



detailed SFHs of galaxies, this section compares these two 
independent methods to estimate the current SFR. 

The purpose of this exercise is three- fold. First, it serves 
as yet another sanity check on the results of the synthe- 
sis. Secondly, it allows us to define in an objective way 
the ratio of current to pa st-average S FR, often referred 
to as Scalo's h parameter (jScalol [l986l ). wh ich is a usefu l 
way to summarize the SFH of galaxies (e.g., 'SandaEe"l986|; 
iBrinchmann et al]|20 04). Finally, defining and calibrating a 
synthesis-based measure of current SFR equivalent to that 
obtained with Ha, allows one to estimate the current SFR 
in galaxies where Ha is not powered exclusively by young 
stars. This turns out to be very useful in studies of AGN 
hosts (|Torres-Papaaui et al.ll2007l . in prep.). 



5.1 Current SFR from Ha luminosity 

For a SFR which is constant over times-scales of the order 
of the lifetime of massive ionizing stars {tion ~ 10 Myr) , the 
rate of H-ionizing photons converges to 



Qh = SFRMh{IMF,Z, 



(8) 



where Nh is the number of hv > 13.6 eV photons produced 
by a SSP of unit mass over its life (in practice, over 95% 
of the ionizing radiation is produced in the first 10 Myr 
of evolution). We computed TVh by integrating the Qnit) 
curves for SSPs using the tables provided by BC03, ob- 
taining Nh = 9.12, 7.08, 6.17, 5.62, 4.47 and 3.16 x 10*^° 
photons Mq^ for Z^, = 0.005, 0.0 2, 0.2, 0.4, 1 and 2.5 Zq, 
respectively, for a lChabrieij (120031 ) IMF between 0.1 and 100 

One in every 2.226 ionizing photons results in emission 
of an Ha photon, almost inde pendently of nebular conditions 
(lOsterbrock fc Ferland|[2006l ). This assumes Case B recom- 
bination and that no ionizing photon escapes the HII region 
nor is absorbed by dust. Adopting the Chabrier IMF and 
the Zq value of A/fr leads to: 



SFRif„ 



2.226Lhc, 



2M0j/r 



-1 / Lhc 
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(9) 



This calibration is strongly dependent on the assumed 
IMF and upper stellar mass limit. Given its reliance on the 
most massive stars, which comprise a tiny fraction of all 
the stars formed in a galaxy, SFRhq involves a large IMF- 
dependent extrapolation, and thus should be used with care. 



5 THE CURRENT SFR 

The most widely employed method to meas ure the "current" 
SFR is by means of the Ha luminosity (|Kennicuttl Il983l ■ 
ll998l : iHopkins et al] |2003| ). We have just devised ways of 
measuring the time dependent SFR which rely exclusively 
on the stellar light, from which one can define a current 
SFR averaging over a suitably defined time interval. Before 
proceeding to the application of these tools to study the 

exclusively with and its associated mass-fraction vector (fl). 
When computing SFRs, however, this di fference must be taken 
into account. Prom the BC03 models for a lChabriei] (|2003|) IMF, 
a 10^'^ yr old population has ~ 2M*, i.e., only half of its 
initial mass remains inside stars nowadays. 



5.2 Current SFR from the spectral synthesis 

The SFR from spectral synthesis is based on all stars that 
contribute to the visible light, and thus should be more rep- 
resentative of the true SFR. We define a mean "current" 
SFR from our time dependent SFHs using equation ((7]l to 
compute the mass converted into stars in the last r years, 
such that 

SFR:(r) = AC-!— (10) 



* Mh is 1.66 times smaller for a Salpeter IMF within the same 
mass limits. 
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is the mean SFR over this period. Because of the discrete 
nature of our base, the function SFR*(r) has a "saw-tooth" 
appearance, jumping every time r crosses one of the t*,j's 
bin borders. 

For the reasons discussed in Section [4.11 it is desirable 
to include components spanning ~ 1 dex in age to obtain 
robust results. Since our base starts at 1 Myr, r ~ 10 Myr 
would be a reasonable choice. This coincides with the mini- 
mum time-scale to obtain SFR* (r) estimates comparable to 
those derived from Lho, which are built upon the assump- 
tion of constant SFR over r > tion ~ 10 Myr. Our base ages 
in this range are = 10, 14, 25, 40 and 55 Myr. 

5.3 Synthesis versus Ha-based current SFRs 

To compare the SFRs given by equations ^ and (llOp we 
must first choose a specific value for r. We do this by cor- 
relating the SFR per unit area obtained with these two es- 
timators, and seeking the value of r which yields the best 
correlation. Surface densities were used to remove the d? fac- 
tors common to both SFRs, thus avoiding distance-induced 
correlations. Data for the SF'"' sample was used in this cal- 
ibration. Also, since Lho refers to the emission from within 
the 3" SDSS fibers, was not extrapolated to the whole 
galaxy in this comparison. 

Fig. [6] shows the results of this exercise. Panel a shows 
the run of the Spearman coefficient {Rs) for different val- 
ues of r, with the best value indicated by an arrow. Given 
the discreteness of our base, any value in the range of the 
ti, — 25 Myr bin yield identically strong correlations (i.e., 
same Rs)- We chose r — 24.5 Myr because this value yields 
zero offset between these two SFRs. This is not a critical 
choice, as values in the whole 10 to 100 Myr range yield 
correlations of similar strength (Fig. |6^). The correspond- 
ing correlations between the synthesis and Ha-based SFRs 
are shown in panels b and c in terms of SFR surface den- 
sities and absolute SFRs, respectively. Robust fits to these 
relations yield slopes very close to unity (1.09 in Fig.[6jD and 
0.94 in Fig.int). 

The rms difference between these two SFR estimators is 
0.3 dex, corresponding to a factor of 2. We consider this an 
excellent agreement, given that these estimators are based 
on entirely different premises and independent data, and 
taking into account the uncertainties inherent to both esti- 
mators. It is also reassuring that r turns out to be compara- 
ble to tion ~ 10 Myr, which is (by construction) the smallest 
time-scale for SFR^fa to be meaningful. That the scatter be- 
tween SFRffa and SFR* is typically just a factor of two can 
be attributed to the fact that we are dealing with integrated 
galaxy data, thus averaging over SF regions of different ages 
and emulating a globally constant SFR, which works in the 
direction of compatibilizing the hypotheses underlying equa- 
tions (O and (fTU)) . 

With these results, we define the ratio of "current" to 
mean past SFR as 

^ _ SFR:(r = 24.5Myr) 
SFR*(r = rG) 

where ra is the age of the oldest stars in a galaxy. In prac- 
tice, since the overwhelming majority of galaxies contain 
components as old as our base allows, the denominator is 
simply M* divided by the age of the Universe, such that 
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Figure 7. Left: Distributions of star formation histories for ^nefc- 
bins A (black) and F (red), as defined in Fig [2] For each bin, we 
show the mean SSFR (soHd fine), the median (dashed) and the 
5 and 95 percentiles of the distributions (dotted). Right: Normal- 
ized distributions of log SSFR for = 25 Myr. Gaussians are 
superimposed to illustrate that the distributions are log normal. 



b is ultimately a measure of the current specific SFR. This 
definition was used in Section [3] where its was shown that b 
decreases by an order of magnitude in the median from the 
top to the bottom of the SF-wing. 



6 THE STAR FORMATION HISTORIES OF SF 
GALAXIES 

Spectral synthesis methods such as the one employed in this 
work have historically been seen with a good deal of skep- 
ticism, best epitomized by ISearld | 1986I ). who, when talk- 
ing about the spectral synthesis of integrated stellar pop- 
ulations, said that "too much has been claimed, and too 
few have been persuaded". Persuading the reader that one 
can nowadays recover a decent sketch of the time-dependent 
SFR in a galaxy from its spectrum requires convincing re- 
sults. In this section we apply the new tools to describe 
SFHs presented above (equations [5] to Q to SF-galaxies in 
the SDSS. As shown below, the temporal-dimension leads 
to a more detailed view of SF-galaxies than that obtained 
with mean ages or current SFR estimates. 



6.1 Distributions of Star Formation Histories 

Our general strategy to explore the statistics of the sample is 
to group galaxies according to certain similarity criteria and 
derive mean SFHs for each group. Since all the results pre- 
sented from Section [6.21 onwards are based on mean SFHs, 
it is fit to first ask how representative such means are of the 
whole distribution of SFHs. 

This is done in Fig [T] where we show the full t*-by-t* 
distribution of SSFR(t*), computed with equation (|6]), for 
two of the six bins in Z„eb defined in Fig [5] bins A and 
F, plotted in black and red, and centered at Z„Eb = 0.31 
and 1.22, respectively. Solid lines indicate the mean SSFR, 
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Figure 6. (a) The solid line shows the Spearman coefficient {Rs) of the E5i?fl(synthesis) X T,sFR(Ha) correlation for different values 
of T in equation JTOj. The dotted line indicates the strength of the SFR{synthesis) X SFR(ifQ:) correlation, (b) Correlation between the 
SFR per unit area obtained through Ha and our synthesis (for r = 24.5 Myr). Units are Mq yr~^ kpc~^ for both axis. The dotted line 
marks the identity line, (c) Correlation between the SFRs derived from equations JTOj and l[9]l. Dashed lines indicate the y{x) and x{y) 
linear regressions, while the solid line shows the bisector fit. 



dashed lines show the median and dotted lines the corre- 
sponding 5 and 95 percentiles of the distributions. The first 
thing one notices in this plot is that the distributions are 
very wide. For instance, for most of the < 1 Gyr range, 
their 5 to 95 percentile ranges span over 2 orders of mag- 
nitude in SSFR. As discussed further below, this is in part 
due to the choice of grouping galaxies by Zneb- Grouping 
by properties more directly related to the SFHs should lead 
to narrower distributions. However, one must realize that 
since galaxy evolution depends on many factors, grouping 
objects according to any single property will never produce 
truly narrow SFH distributions. Secondly, the distribution 
of SSFR values at any is asymmetric, as can be seen by 
the fact that the mean and median curves differ. In fact, as 
illustrated by the right panel in Fig [71 these distributions 
are approximately log-normal, indicating that SFHs result 
from the multiplication of several independent factors, as 
qualitatively expected on physical grounds. 

Despite their significant breadth and overlap, it is clear 
that the SSFR-distributions for Z^eb bins A and F in Fig [7] 
are very different, particularly at low f*. This is confirmed 
by KS tests, which show that these two distributions are 
undoubtedly different. In fact, for any pair of Znet bins the 
distributions differ with > 99% confidence for any , 

In what follows we will present only mean SFHs, ob- 
tained grouping galaxies according to a subset of the avail- 
able physical parameters. Whilst their is clearly more to be 
learned from the SFH distributions discussed above, this is 
a useful first approach to explore the intricate relations be- 
tween galaxy properties and their SFHs. 

6.2 Trends along the SF-wing 
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We start our statistical study of galaxy SFHs grouping 
galaxies in the six Zn^i, bins defined in Section [21 As shown 
in Fig[2l Z^et traces the location of a galaxy along the SF- 
wing in the BPT diagram. 

Fig [SI shows the derived SFHs for the A-F bins in 
four different representations, from top to bottom: Xs , SFR, 
SSFR and 7)J as a function of stellar age t*. Each line rep- 



Figure 8. Mean star formation histories for the 6 different 
•^nefc-bins defined in Fig [2] in four different representations: (a) 
smoothed population vector, x7(tj,), (b) SFR(f*), (c) SSFR(ty,), 
and (d) 
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resents a f*-by-t* average over all galaxies in the bin. The 
plots show that young populations are present in a propor- 
tion which increases systematically as Zneb decreases. This 
is evident, for instance, in the a;s(t*) panel, which shows 
how Z„ei,-related age distributions combine to produce the 
correlation between (logt+)i and Z^sb depicted in Fig. [Sf. 

The SFR(t*) curves (panel b) show that SF-galaxies of 
different Zn^b differ more in their past SFR, low Zneb having 
SFRs about 100 times lower than those of high Zneb a few 
Gyr ago. In the more recent past, all curves converge to 
SFRs of a few MQ/yr. At first sight, this convergence seems 
to be at odds with the fact that galaxies with low and high 
Zneb differ by about one order of magnitude in the median 
Lhc (Table [!}, and thus should differ by a similar factor in 
the recent SFR. In fact, there is no contradiction, since what 
needs to be considered when comparing the synthesis-based 
SFR with that derived from Ha is the mean SFR over scales 
of at least 10 Myr, and these are clearly smaller for low Zneb 
galaxies than for those of higher Zneb- As shown in Fig [B] 
Ifa and synthesis based SFRs agree very well. Ultimately, 
the apparent coincidence of mean SFR curves of different 
Zneb is due to the fact that the relation between recent SFR 
and Zneb is a relatively weak and scattered one (Fig. [SJi) , 
such that along the whole SF-wing one may find galaxies 
that transform a similar amount of gas into stars per year. 

The clearest separation between SFHs of galaxies of dif- 
ferent Zneb is provided by the SSFR(t*) curves. At ages > 
a few Gyr all SSFR curves merge. This behavior is a conse- 
quence of the fact that most of the stellar mass is assembled 
early on in a galaxy's history, irrespective of Zneb or other 
current properties. With our Alogt* = 1 dex smoothing, 
this initial phase, over which J SFRdf* ~ M^, becomes a 
single resolution element in the SFR curves. Division by 
to produce a specific SFR (equation[6]) then makes all curves 
coincide. At later times (smaller t*), however, the curves di- 
verge markedly, with the lowest and highest Zneb groups 
differing in SSFRs by ~ 2 orders of magnitude nowadays. 
This confirms that the relation between recent and past star- 
formation is a key-factor in distributing galaxies along the 
SF-wing in the BPT diagram (Fig. [5^). 

Yet another way to visualize the SFH is through the 
mass-assembly function defined in equation ((Tjl. Though 
ViKti,) is computed with the raw (unsmoothed) population 
vector, for presentation purposes we apply a FWHM = 0.2 
dex gaussian in logt*, just enough to smooth discontinuities 
associated with the discrete set of t«,j's in our base. Fig|8]l 
shows the results. This is essentially a cumulative represen- 
tation of the same results reported in Fi^H};, namely, that 
low Zneb galaxics are slower in assembling their stars. This 
plot is however better than the previous ones in showing that 
despite these differences, all galaxies have built up most of 
their stellar mass by = 1 Gyr. 

These encouraging results indicate that synthesis meth- 
ods have evolved to a point where one can use them in con- 
junction with the fabulous data sets currently available to 
sketch a fairly detailed semi-empirical scenario for galaxy 
evolution. In the next section we walk a few more steps 
in this direction by inspecting how astrophysically plausi- 
ble drivers of galaxy evolution relate to the SFHs recovered 
from the data. 



6.3 Star Formation Histories and Chemical 

Evolution in Mass and Surface-Density bins 

The value of grouping galaxies by Zneb is that it maps SFHs 
to a widely employed diagnostic tool: the BPT diagram. 
Yet, present day nebular abundance is not a cause, but a 
consequence of galaxy evolution. In this section we leave 
aside our focus on the BPT diagram and group galaxies ac- 
cording to properties more directly associated to physical 
drivers of galaxy evolution. Two natural candidates are the 
mass {Mi,) and surface mass-density (E*). Like Zneb, both 
Mi, and can be considered the end product of a SFH, yet 
they are clearly more direct tracers of depth of the potential 
well and degree of gas compression, two key parameters af- 
fecting physical mechanisms which regulate galaxy evolution 
(Schmidt., 195a : ,Tinslev,198a : ^Kennicutt 1998 ). 

Fig |9] shows our different representations of the SFH 
of SF galaxies for five 1 dex-wide mass bins centered at 
log A/* /A/0 = 7.5 • • • 11.5. Given that Af« and Zneb are re- 
lated, the overall evolutionary picture emerging from this 
plot is similar to the one obtained binning galaxies in Zneb, 
i.e., massive galaxies assemble their stars faster than low- 
mass galaxies. On the whole, Fig [9] provides a compelling 
visualization of galaxy downsizing. 

The most noticeable difference with respect to Zneb- 
binned results is on the absolute SFR curves (compare Figs 
(Sja and[9l3). This difference is rooted in the fact that galax- 
ies of similar Zneb span a much wider range of SFRs than 
galaxies of similar Mi,. This can be illustrated focusing on 
recent times, and inspecting the Lna-Zneb relation (Fig. 
[SJi), with the understanding that Lhc can be read as the 
current SFR (Fig. [6]). Despite the statistically strong cor- 
relation {Rs — 0.49), the typical 5-95 percentile range in 
Lho for a given Zneb is ~ 2.1 dex, comparable to the full 
dynamic range spanned by the data (2.4 dex over the same 
percentile range). In other words, the relation has a large 
scatter, and hence Z^eb-binning mixes objects with widely 
different SFRs, explaining why the SFR(f*) curves in FiglSja 
tend to merge at low t*. The I/Ha-Af* relation (not shown), 
on the other hand, is stronger {Rs — 0.68), partly due to 
the factors in common to absolute SFRs and Af*. Group- 
ing by Mi, then selects galaxies in narrower SFR ranges, 
producing the well separated curves seen in Fig[U]D- 

Results grouping galaxies according to E* are presented 
in Fig [lOl Since E* and A/* correlate very strongly {Rs — 
0.73 in our sample), the results are similar to those obtained 
gr oupin g galaxies by their stellar mass. iKauffma nn et al.l 
(|2003bi ). based on an analysis of two SFH-sensitive spec- 
troscopic indices (namely _D„(4000) and H5a), propose that 
E* is more directly connected to SFHs than Af* . This is not 
obviously so comparing Figs [9] and [10] A more detailed, mul- 
tivariate analysis is needed to evaluate which is the primary 
driver of SFHs. 



7 SELECTION EFFECTS AND MODELLING 
CAVEATS 

This section deals with the effects of sample selection, syn- 
thesis ingredients and model assumptions on our results. 
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Figure 9. As Fig |8] but binning SF galaxies by their stellar 
mass, using five 1 dex wide bins, centered at (from bottom to 
top in panel b) logAf*/M0 = 7.5, 8.5, 9.5, 10.5 and 11.5, which 
contain 542, 4057, 26700, 47153, 3808 galaxies, respectively. 
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Figure 10. As Fig [S] but binning SF galaxies by their stellar 
surface densities, with five 0.5 dex wide bins centered at log S* = 
7.25, 7.75, 8.25, 8.75 and 9.25 Af© kpc-^, containing 1477, 12554, 
37177, 27742, 3046 galaxies respectively. 



T.l Selection effects 

We now study to which extent the mean SFHs derived in 
the last section are affected by the way we have defined SF 
galaxies. We address this issue recomputing mean SFHs for 
samples constructed with alternative selection criteria, and 
comparing to the results obtained with our default sample. 

We first ask how our emission line and continuum S/N 
cuts infiuence our results. Figs lllb . and b show the average 
mass assembly and SSFR functions for the high-quality SF*^' 
sub-sample defined in Section [?!3l The results for this better- 
data sub-sample are very similar to those obtained with the 
full SF sample. The SSFR curves in recent times are skewed 
to slightly higher rates, which reflects the fact that objects 
in the SF'"' sample are slightly younger than those in the 
full SF sample, as shown in Fig|3l 

Our BPT-based selection of SF galaxies used the di- 
viding line proposed in SEAGal III, whic h is more restric- 
tive th an the empirical line proposed by iKauffmann et ahl 
lj2003d '). We define the SF*'' sam ple as the 111026 galaxies 
classified as SF according to the IKauffmann et al.1 (|2003cl ') 
line. Fig. Illb and d show that the SFHs for Zne.b bins in 
this sample are nearly indistinguishable from those obtained 
with the SEAGal classification criterion. 



Another concern is the inclination effect. The spectra of 
edge-on objects are biased by the metal poorer outer parts 
of the galaxies, leading to an underestimation of Zneb- This 
may lead us to place a galaxy in a lower Z^eb bin than it 
would if it were seen face on, possibly affecting the mean 
SFH in that bin. To investigate this effect we have defined a 
sub-sample of nearly face-on galaxies, SF^° (6842 objects), 
selecting by the inclination parameter, h/a > 0.9. Figs lllfe 
and f show that the SFHs derived with this sample are prac- 
tically the same as for the full sample. 

Aperture effects are a co mmon source of co ncern in 
studies of SDSS spectra (e.g., iGomez et a"l]|2003l : SEAGal 
I). To investigate how such effects impact upon our SFHs 
we defined two samples: the SF^ sample, which comprises 
58153 SF galaxies with z > 0.05 (as opposed to 0.002 for 
the full SF sample), and the SF"*" sample, comprising only 
the 1096 objects with more than half of their z band lumi- 
nosity inside the fiber. Both criteria exclude proportionately 
more the population of low Z^eb, low Af^. galaxies (distant 
galaxies of this category are not present in the SDSS, due to 
limiting magnitude) . Accordingly, the change in SFHs is only 
noticeable for the lowest Z^eb bins, as shown in Figs lllb -i. 
In particular, Fig llll shows that the mean SFH for bin C in 
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Figure 11. Average mass assembly histories (77^, top panels) and specific star formation rate (SSFR, bottom) histories for Z„e(j-bins B 
to F, color-coded as in Fig [2] Solid lines show the curves for different sample definitions. For comparison, the dotted lines in all panels 
show the curves for bins B (lower curves; magenta) and F (upper curves; red) of the full SF sample. 



the SF"^ sample matches that of bin B in the fuU sample. 
Of all selection-induced changes discussed here, this is the 
largest one; yet, all it does is to shift the SFHs from one 
group of galaxies to the adjoining one. 

To summarize, selection criteria may influence the de- 
rived mean SFHs in quantitative terms, but do not modify 
the relative pattern of mean SFHs of galaxies in different 
Z„ei) bins. The same applies to grouping galaxies accord- 
ing to properties other than Z^eb- The general trends in the 
SFH as a function of global galaxy properties obtained in 
this work are therefore robust against variations in the sam- 
ple selection criteria. 

7.2 Experiments with different models 

One should also ask to which extent our results are robust 
against changes in the base of evolutionary synthesis models. 
While answering this question requires an in-depth study far 
beyond the scope of this paper, we believe this has a much 
larger impact on SFHs than selection effects. 

Panter et al. (2007) reported results of MOPED ex- 
periments usin g spec tr al models from different s ources 
ll Jimenez et all 20041; iFioc fc Rocca-Volmerangd 1 19971 ; 
iMarastonI |2005| ; Isruzual A. fc Chariot! Il993l ; and BC03), 
all sampled at AA = 20 A. For the 767 galaxies in their 
randomly selected test sample, the resulting mean star- 
formation fractions (analogous to our /2 vector) differ by 
factors of a few for the youngest and oldest ages, and close 
to a full order of magnitude for ^ 0.1-1 Gyr. Recovering 
SFHs in this intermediate a ge ran ge is particularly hard, as 
discussed bv iMa this et al.l (|2006l ). Indeed, the experiments 
reported by Panter et al. often find a suspiciously large 1 
Gyr component. 

The behaviour of our starlight fits is also anomalous 
in this age range. This is clearly seen in the mean SFHs 
shown in Figs IHl llOl particularly in the Xs representation. 



which shows a hump at ~ 1 Gyr. Interestingly, starlight 
experiments with a base of evolutionary synthesis models 
using the MILES library of ISanchez-Blazguez et al] (|2006l . 
instead of the STELIB library used in the base adopted for 
this and previous SEAGal studies) do not produce this hump 
at ~ 1 Gyr. In fact, the whole mean SFHs derived with this 
new set of models differs systematically from those shown in 
Figsl Mini The mass assembly function rjt{tt), for instance, 
rises more slowly and converges at somewhat smaller than 
those obtained in this work. There are also systematic dif- 



ferences in stellar extinction, which comes out AAt 



0.4 



mag larger with the MILES models. Extensive tests with 
these new bases, including different prescriptions of stellar 
evolution as well as different spectral libraries, are underway, 
but these first results show that significant changes can be 
expected. 

Reassuringly, however, these same experiments show 
that the pattern of SFHs as a function of Zmt, Mi, and 
E* reported in this paper does not change with these new 
models. In any case, these initial results provide an eloquent 
reminder of how dependent semi-empirical SFH studies are 
on the ingredients used in the fits. 

T.3 Experiments with differential extinction 



We saw in Section 12.5.21 that the nebular and stellar exct- 
intions are strongly correlated, but with Ay"^ ~ twice Ay. 
This indicates that the uniform stellar extinction used in our 
fits is not adequate to model star-forming regions, which 
should be subjected to a similar extinction than the line- 
emiting gas, ie., j45>(t* < lO^yr) ~ Ay"''. It is therefore fit to 
ask whether and how such a differential extinction affects 
our general results. 

Given the difficulties in recovering reliable population 
dependent Ay's from spectral synthesis in the optical range 
alone, we address this question by postulating that popula- 
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Figure 12. Comparison between mean star formation histories 
for the SF'"' sample modelled with only one stellar extinction 
(dotted lines) and with differential extinction (solid lines). For 
clarity, only Z„ei,-bins B, D and F are drawn, and color-coding is 
the same as in Fig [2] Each panel shows a different representation 
of SFHs: (a) smoothed population vector, x7(t*), (b) SFR(t*), 
(c) SSFR(t*), and (d) r^(t*). 



tions younger than 10^ yr are extincted by 0.34 + 2.2SAy, 
i.e., the empirical Ay^^{Ay) relation found in Section P.5.21 
with Ay now denoting the extinction to > 10^ yr starsjj The 
17142 galaxies in the SF'"' sub-sample were re-fitted with 
this more realistic modified recipe for extinction effects. 

Qualitatively, one expects that forcing a uniform Ay fit 
to a galaxy where the young stars suffers more extinction 
than the others should lead to an overestimation of the age 
of the young population. This older and thus redder young 
population would compensate the mismatch in Ay. Allowing 
for Aylt < lO'^ yr) larger than the Ay of the t > 10^ yr 
populations makes it possible to fit the same spectrum with 
younger and dustier populations. Hence, the recent SFRs 
should increase. 

This expectation is fully confirmed by these new fits. 
Fig. [12] shows a comparison between three ^„e6 bins (B, D 
and F) for the old and new fits of the SF'"' sample. One sees 
that the new average SFR and SSFR curves are shifted by ~ 
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0.3 dex upwards in the < 10^ yr range with respect to the 
ones obtained with a single extinction. The rearrangements 
in the population vector tend to be in the sense of shifting 
some light from old populations to the < 10^ yr components. 

Not surprisingly, the properties which change most are 
those directly related to the strength of the young popu- 
lation, such as current SFR, which increases by 0.3 dex 
on average, and the mean stellar age, which decreases by 
A(logf*)L ~ 0.1 dex. The changes in other global proper- 
ties such as Ay and A/* are much smaller than this. 

These experiments are obviously a simplification of the 
problem of dust distribution in galaxies, yet they suggest 
that the choice of the extinction modelling can have non- 
negligible effects on the derived SFH curves. On the whole, 
however, the qualitative pattern of SFHs as a function of 
Z„eb or other variables stays the same. As found in the sam- 
ple selection studies, quantitative changes are at best equiv- 
alent to moving from one bin to the next, so that our general 
conclusion does not change. 



8 SUMMARY 

In this paper we have studied physical properties of 82302 
normal star-forming galaxies from the SDSS DR5, by us- 
ing results from our stellar population synthesis code, 
STARLIGHT, and our emission-line measurement algorithm. 

Before reviewing our main results, we highlight some 
aspects of this study which have relatively little impact 
upon our general conclusions, but represent significant re- 
finements in our methodology. 

(i) We have detected a systematically overestimated con- 
tinuum level around H/3, whose origin we tentatively at- 
tribute to deficiencies in the STELIB calibration in this 
range. Gaussian fits to the H/3 emission which disregard this 
offset tend to underestimate the line flux by 4% on average, 
and ~ 7% in the case of weaker lines. These are relatively 
small, yet systematic effects, which propagate to estimates 
of nebular extinction, metallicities and galaxy classification. 

(ii) SF galaxies were selected according to the theoretical 
criterion proposed in SEAGal HI, which minimizes contam- 
ination by AGN emission. 

(iii) Nebular extinctions and metallicities were derived 
self-consistently, allowing for the metallicity dependence of 
the Balmer decrement. Five different reddening laws were 
explored, but found to produce equally good spectral fits 
and relatively small differences in derived physical proper- 
ties. 

(iv) We have confirmed the strong correlation between 
Ay"'' and Ay found in SEAGal I. We have also identified a 
strong correlation between the strength of the ISM compo- 
nent of the NaD absorption doublet and the amount of dust 
derived from the synthesis. 

(v) Different recipes for nebular metallicity estimates 
were tried. Some of them proved not to be adequate for 
this study, either because of the lack of spectral data (e.g., 
measures of [ArIII]A7135 and [OIIIJ4363 emission lines were 
available for few objects), or because such calibrations were 
only valid in the low-Z„eb regime, thus encompassing a 
very small fraction of objects from our sample. Therefore, 
throughout our analy s is we use the O3N2 index and the cal- 
ibration bv lStasihskal l|2006l ) to measure the nebular metal- 
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licity. Although this is not a rehable caUbrator at the lowest 
metallicities, it is good enough for our analysis in Znet-bins. 
Furthermore, it has the nice virtue of being directly related 
with the position of the objects in the BPT diagram. 

We now summarize results related to the main goal of 
this paper, namely, to investigate the SFH of galaxies along 
the SF wing in the BPT diagram. In practice, this means 
studying how SFHs change as a function of nebular metal- 
licity, even though Z^eb is more a product than a cause of 
galaxy evolution. 

(i) We started our study with a traditional analysis, cor- 
relating Zneb with several physical and observed properties. 
This analysis confirms results obtained directly or indirectly 
in the past by other works, such as relations between the 
nebular metallicity and galaxy luminosity, mass, dust con- 
tent, mean stellar metallicity and mean stellar age. 

(ii) Formalism towards a time-dependent analysis was 
then presented. Simple ways to compress the output of our 
stellar population synthesis code were proposed. These are 
based either on a posteriori smoothing of the age distri- 
bution, which allows the derivation of time dependent star 
formation rates, or a cumulative mass assembly history. 

(iii) As a first application of this time dependent descrip- 
tion of SFHs we computed the current SFR obtained from 
our spectral fits. The resulting values of SFR* agree very well 
with more traditional estimates based on the luminosity of 
Hq. The scatter between SFR* and SFRh^ is just a factor of 
2, despite the differences in the underlying assumptions and 
sensitivity to the IMF. This result strengthens confidence in 
our method, and, more importantly, opens the possibility of 
measuring current SFRs in galaxies hosting AGN, where Ha 
method does not apply. 

(iv) Fully time dependent SFHs were then derived group- 
ing galaxies into six Z^eb bins spanning the entire SF wing 
of the BPT diagram. Mean SFHs for each of these bins 
were presented in four different representations; (a) the 
smoothed population vector, (b) the star formation 
rates SFR(tt), (c) specific star formation rates SSFR(t*), 
and (d) mass-assembly histories , rj1{ti,). 

(v) We found that SFHs vary systematically along the SF 
sequence. Though all galaxies assembled the bulk of their 
stellar mass over 1 Gyr ago, low Zneb systems evolve at a 
slower pace. Galaxies at the tip of the SF wing have current 
specific SFRs about 2 orders of magnitude larger than the 
metal rich galaxies at the the bottom of the BPT diagram. 

(vi) At any given time, the distribution of SSFRs for 
galaxies within a Z,jei,-bin is quite broad and approximately 
log-normal. 

(vii) We performed the same SFH study grouping galax- 
ies by their stellar mass and surface mass density. Given the 
existence of Z„ei,-M*-E* relations, the overall picture is ob- 
tained grouping by Z^eb- Thus, low M* (low E*) systems 
are the ones which evolve slower, with current SSFRs much 
larger than more massive (dense) galaxies. 

(viii) Finally, we have analysed a number of selection and 
modelling effects that might bias our results, and show that 
while they may affect the derived SFHs quantitatively, the 
organization of SFHs as a function of Zneb, Mi,, re- 
mains the same. Experiments with new evolutionary syn- 
thesis models and differential extinction fits were reported 



and found to lead to substantially different SFHs, yet pre- 
serving this same overall pattern. 
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